Comment on "Quantum Phase Transition of the Randomly Diluted 
Heisenberg Antiferromagnet on a Square Lattice" 
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In H , Kato et al. presented quantum Monte Carlo re- 
sults indicating that the critical concentration of random 
non-magnetic sites in the two-dimensional antiferromag- 
netic Heisenberg model equals the classical percolation 
density; p c — p* « 0.407254. The data also suggested 
a surprising dependence of the critical exponents on the 
spin S of the magnetic sites, with a gradual approach to 
the classical percolation exponents as S — > oo. This be- 
havior is contrary to well established notions of universal- 
ity, according to which the exponents should depend only 
on symmetries and dimensionality. I here argue that the 
exponents in fact are ^-independent and equal to those 
of classical percolation. The apparent S'-dependent be- 
havior found in |ffl] is due to temperature effects in the 
simulations as well as a quantum effect that masks the 
true asymptotic scaling behavior for small lattices. 

I discuss simulation results for the S = 1/2 case ob- 
tained using the stochastic series expansion method on 
LxL lattices at the percolation density. Fig. [l]a shows re- 
sults for the disorder-averaged staggered structure factor 
(the equal-time spin correlation function at wave-vector 
q = (7r,7r), defined as in [jjj) versus inverse temperature 
for three different system sizes (J is the Heisenberg cou- 
pling). The same configurations of magnetic sites were 
used at all temperatures and the statistical errors are 
therefore correlated and completely dominated by fluctu- 
ations between different site configurations (several hun- 
dred random site configurations were used for each L). 
In lattices with L < 48 were studied at tempera- 
tures down to J/1000. Fig. |l|a reveals that for L — 48 
the structure factor has not yet converged to its asymp- 
totic T = value at this temperature. A careful sta- 
tistical analysis, taking into account the covariance be- 
tween data at different temperatures, gives that S(ir, tt) 
at J/T = 1024 is (95.0 ± 0.5)% of the almost saturated 
value at J/T = 8192. It is thus likely that the finite-size 
scaling in ffl is to some extent affected by temperature. 

Temperature effects alone cannot explain the anoma- 
lous behavior of the critical exponents found for small S 
in Q . In order to better understand the scaling behavior 
of S(tt, 7t) it is useful to also study the cluster structure 
factor 5 c i(7r, tt) defined on the largest connected cluster of 
magnetic sites. The corresponding cluster magnetization 
is given by = (3S c i(tt,tt)/N c i), where N c \ is the size 
of the cluster and ( ) denotes the average over random 
site configurations. In the thermodynamic limit the clus- 
ters have fractal dimension D = 91/48 and in [|] it was 
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FIG. 1. (a) Temperature dependence of the structure factor 
for different system sizes, (b) Size-dependence of the cluster 
magnetization for L in the range 8,. ..,48. 


argued that the spin correlations on these clusters de- 
cay algebraically with an ^-dependent exponent. Fig. ^> 
shows m^j versus 1/L D ^ 2 along with a quadratic fit. The 
cluster magnetization clearly scales to a finite value as 
L — > oo and therefore, contrary to the claim of [0, the 
fractal clusters for S = 1/2 are long-range ordered. The 
scaling of the full S(ir, tt) must then be classical and, as 
was also concluded in |jj , the critical density p c = p*. 

The reason why the classical exponents are not seen 
in the full staggered structure factor considered in Q is 
that the cluster magnetization decays considerably with 
increasing lattice size for the sizes used, as can be seen 
in Fig. |l|b. This partially compensates for the growth of 
the fractal clusters with L, leading to a slower growth of 
S(ir, tt) with L than in a classical system for which 
is constant. The asymptotic behavior can be observed 
only when L is sufficiently large for the relative size cor- 
rections to m^j to be small. This quantum mechanical 
effect makes it extremely difficult to directly observe the 
true critical behavior for small values of S. 
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